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The compact X-ray binary system Cyg X-3 has been recently discovered as a source 
of GeV 7-rays by the AGILE and the Fermi satellites. It shows emission features in 
the GeV 7-rays similar to other 7-ray binaries which were also observed in the TeV 
, 7-rays (LS 5039 and LSI +61 303). The question appears whether Cyg X-3 can be 

also detected in the TeV 7-rays by the Cherenkov telescopes. 

Here we discuss this problem in detail based on the anisotropic inverse Compton 
(IC) e ± pair cascade model successfully applied to TeV 7-ray binaries. We calculate the 
7-ray light curves and 7-ray spectra expected from the cascade process occurring inside 
the Cyg X-3 binary system. It is found that the 7-ray light curves at GeV energies 
can be consistent with the 7-ray light curve observed by the Fermi for reasonable 
parameters of the orbit of the injection source of relativistic electrons. Moreover, we 
show that in such a model the sub- TeV 7-ray emission (above 100 GeV) is expected to 
CNI ' be below sensitivities of the present Cherenkov telescopes assuming that electrons are 

accelerated in Cyg X-3 to TeV energies. The next stage Cherenkov telescopes (MAGIC 
CO , II, HESS II) should have the energy threshold in the range 20-30 GeV, in order to 

have a chance to detect the signal from Cyg X-3. Otherwise, the positive detection of 
7-rays at energies above a few tens of GeV requires a telescope with the sensitivity of 
~ 0.1% of Crab Units. We conclude that detection of sub- TeV 7-rays from Cyg X-3 
CN ■ by on-ground telescopes has to probably wait for the construction of the Cherenkov 

O Telescope Array (CTA). 
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1 INTRODUCTION cm -2 s _1 (Abdo et al. 2009a). Interestingly, the 7-ray light 

curve at GeV energies from Cyg X-3 shows general features 

The enigmatic binary system, Cyg X-3, has been claimed that are <l uite similar to those observed recently from the 

as a powerful 7-ray source at different energies since the TeV ^" ra y binaries LS 5039 and LSI 303 +61, which were 

first satellite and on-ground observations starting from 70- also observed by the Fermi-LAT at GeV energies (Abdo et 

80-ties. The reports on the MeV-GeV 7-ray emission were al - 2009b,c). Therefore, future observations of Cyg X-3 with 

diverse. A marginal detection has been claimed by the SAS-2 the modern Cherenkov telescopes seem to be well motivated, 

and EGRET (Lamb et al. 1977, Mori et al. 1997). How- In fact, during last years the MAGIC telescope has observed 

ever, the negative result has been reported by the COS °yg X " 3 for ~ 70 hours in diff erent emission states. No pos- 

B (Hermsen et al. 1987). Also early optimistic results re- itive si S nal has been reported up to now (Saito et al. 2009). 

ported by the Cherenkov telescopes were not confirmed (e.g. The u PP er limits are on the level of 1% of the Crab Unit 

Weekes 1988, 1992; Chadwick et al. 1990). Only recently, the (OIL) above ~ 250 GeV. 

AGILE and Fermi satellites reported positive transient de- The observed modulation of GeV 7-ray emission from 

tection of Cyg X-3 at GeV energies (Tavani et al. 2009, Abdo Cyg X-3 strongly suggests that photons has to originate in 

et al. 2009a). The 7-ray emission from Cyg X-3 has been dis- the radiation process which occurs inside the binary system, 

covered before the major radio flares (Tavani et al. 2009). The most likely scenario is the interaction of electrons with 

It shows a modulation with the period of the binary system the anisotropic radiation of the WR type companion star. 

(Abdo et al. 2009a). The emission above 100 MeV is well de- Electrons are accelerated at some physical process either 

scribed by a single power law with differential spectral index close to accreting object or in the jet or in the vicinity of a 

— 2.7+0. 05(stat)±0.20(syst) and the peak flux ~ 2xl0 -6 ph pulsar or at the shock wave between the jet and the stellar 
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wind (as postulated by models proposed for TeV 7-ray bi- 
naries). However, due to very strong radiation field created 
by such luminous star inside compact binary Cyg X-3 (sur- 
face temperature T* ~ 10 5 K, stellar radius i?* ~ 2 x 10 11 
cm, and separation of the compact object D ~ 7- 
rays should initiate IC e pair cascades as considered for 
TeV 7-ray binaries (e.g. Bednarek 2000, Orellana et al. 2007, 
Khangulyan et al. 2008). Such IC e* pair cascade scenario 
has been already investigated also for the Cyg X-3 binary 
system (Bednarek 1997 and Sierpowska & Bednarek 2005). 
It has been shown in those papers that the optical depths for 
7-rays created inside the Cyg X-3 binary system are huge 
(they can reach values up to a few hundred in specific di- 
rections). Therefore, escape of 7-rays with energies above a 
few tens of GeV is not very likely. In this paper, we perform 
detailed calculations of the 7-ray spectra escaping towards 
the observer located at different directions in respect to the 
orbital plane of the Cyg X-3 binary system for different loca- 
tions of the source of energetic leptons within the binary. We 
calculate the expected 7-ray light curves at energies suitable 
for the satellite and Cherenkov telescopes. 



2 PHASE DEPENDENT GAMMA-RAY 
OPTICAL DEPTHS 

The optical depths of 7-ray photons in the anisotropic ra- 
diation of the massive star has been calculated in general 
case (arbitrary injection place within the binary) for the 
star with the parameters of the WR star in Cyg X-3 by 
Bednarek (1997, see Fig. 2 in that paper). Here we calculate 
the 7-ray optical depth light curves for the case of 7-rays 
injected from a point source on a circular orbit around the 
star with a radius, Rj. Since the inclination angle of the Cyg 
X-3 binary system is expected in the range i = 30° — 60°, 
we show how the optical depths change with the phase of 
the binary system for the limiting values of the inclination 
angle. We apply the following parameters for the WR star in 
Cyg X-3: its radius R* = 2 x 10 11 cm, surface temperature 
T* = 9 x 10 4 K and the radius of the orbit of the compact 
object a = 2.25-R*. These parameters are consistent with the 
limits put by Cherepashchuk & Moffat (1994). Note that the 
injection place of 7-rays does not need to be identical with 
the location of the compact object in this binary system. 

We consider two general geometrical situations for the 
location of the injection source of 7-rays: 

(i) It is in the plane of the binary system but at different 
distances from the companion star (in the range i? 7 = 1.5 — 
20i?*). Such geometry can give us an idea about the 7-ray 
emission in the case of acceleration regions extended in the 
plane of the binary system (e.g. the shock structure created 
in the collision of the stellar and pulsar winds). 

(ii) It is above the plane of the binary system for the 
distance of the compact object from the companion star (in 
the range D = 1 — 20i?*). Such calculations are suitable in 
the case of the acceleration region extended in the direction 
of the jet which is perpendicular to the plane of the binary. 

The optical depths for 7-rays propagating in the 
anisotropic radiation of the companion star from the injec- 
tion place up to the infinity for these two situations are 
shown in Figs. 1 and 2. In the case of a circular orbit with 



different radii, we note clear variability of the optical depths 
for 7-rays with different energies. The general tendency is 
that the level of variability of the optical depth decreases 
with energy of 7-ray photons in the range of ~33 GeV to 
~1 TeV. The 7-ray optical depths change more drastically 
with the phase of the binary system for larger inclination an- 
gles. As expected, the optical depths reach the lowest values 
when the injection source is in front of the companion star 
(phases close to 0.5). Optical depths can reach values below 
unity at certain region around the phase 0.5. This region be- 
comes less extended for larger energies of 7-ray photons. In 
the case of the orbit with the radius _R 7 = 1.5r*, the optical 
depths decrease for small and large phases due to the eclipse 
of the 7-ray source by the companion star (in this case the 
optical depths are calculated only to the stellar surface). 

The 7-ray optical depth light curves calculated for dif- 
ferent distances above the plane of the binary system show 
also local minima at phases corresponding to the location 
of the injection place in front of the companion star (see 
Fig. 2). These minima become at first deeper with increas- 
ing distance of the injection place from the orbital plane. 
The appearance of these features is related to the decrease 
of the angle between direction towards the observer and the 
direction towards the injection place of 7-rays (counted from 
the stellar center) for larger distance, D, from the orbital 
plane. However, for large distances, D, the dips in the 7- 
ray light curves become shallower and the general level of 
the optical depth becomes lower. As a result, modulation 
of the optical depths with the period of the binary system 
becomes relatively low. Based on these calculations, we con- 
clude that the largest contribution to the modulation of the 
7-ray emission with the period of the binary system comes 
from the intermediate distances from the orbital plane. How- 
ever, the distance, from which the maximum contribution to 
the modulation of the 7-ray signal is seen, depends also on 
the inclination angle of the binary system. 

Based on the calculations of the phase dependent opti- 
cal depths for 7-ray photons with different energies, we can 
conclude that the cascading effects should be very important 
in the case of Cyg X-3 binary system for both considered in- 
clination angles. However, the 7-ray spectra escaping to the 
observer have a chance to extend up to the sub- TeV energies 
only for a range of phases around tp — 0.5, corresponding 
to the location of the injection region in front of the com- 
panion star. On the other hand, when the injection source 
is behind the companion star, the optical depths are clearly 
above unity for most of the considered locations of the injec- 
tion source and energies of 7-ray photons. For these phases, 
7-ray emission should be limited only to energies below a 
few tens of GeV. In the next section we discuss the 7-ray 
emission features in terms of a simple scenario in which rel- 
ativists electrons are injected isotropically from a point like 
source located within the binary system. 



3 GENERAL SCENARIO FOR GAMMA-RAY 
PRODUCTION 

The observed modulation of GeV- TeV 7-ray emission with 
the period of the binary systems (LS 5039, LSI 303 +61, 
Cyg X-3) indicates that the production of 7-rays occurs in- 
side the volume of the binary system. The companion stars 
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Figure 1. The optical depths for 7-ray photons injected at a specific place inside the binary system for energies Ej = 33 GeV (dashed 
curves), 100 GeV (solid), 330 GeV (dot-dashed), and 1 TeV (dotted) as a function of the phase of the binary system Cyg X-3. The phase 
zero corresponds to the location of the injection source behind the massive star. Figures on the left are obtained for the inclination of the 
binary system equal to i = 30° and on the right for i = 60°. The injection source is on the circular orbit within the plane of the binary 
system but at different distances from the companion star: i? 7 = 1.5R* (a and f), 2.25R* (t>,g), 5-R* (c,h), 10R* (d,i),20-R* (e,j). Thin 
dotted line marks the optical depth equal to unity. The companion star has the radius fl t = 2x 10 11 cm and the surface temperature 
T* = 9 X 10 4 K (as expected for the companion star in the Cyg X-3 binary system). 



in these binaries are very luminous which again suggest that 
emission process is related to strong radiation field created 
by these stars. At the conditions envisaged above, the pro- 
duction of 7-rays likely occurs in the cascade process since 
the cross sections for the inverse Compton process and the 



e pair creation in collisions of 7-rays with the stellar low 
energy photons are comparable. In fact, relativistic hadrons 
might also contribute to the 7-ray flux provided that they 
find enough target provided by the dense stellar winds or 
the accretion disk around a compact object. Due to the ob- 
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Figure 2. As in Fig. 1 but for the case of the 7-ray injection source located in the jet launched from the compact object perpendicularly 
to the plane of the binary system. The compact object is at the distance of a = 2.25R* from the companion star. Specific figures show 
the optical depths for distances along the jet equal to: D = li?* (a,f), 3R* (t>,g), 5R* (c,h), 10R* (d,i), 20.R* (e,j). 



served modulation of the 7-ray signal, the radiation pro- 
cesses likely occur within the binary system. The specific 
radiation mechanism which is consistent with the observed 
features of the 7-ray emission from the binary systems has 
been already studied before the discovery of TeV 7-ray emis- 
sion from massive binary systems (see e.g. Bednarek 2000, 
Sierpowska & Bednarek 2005). It assumes that relativistic 
electrons are accelerated inside the binary system in one of 
the popular general picture, i.e either by the pulsar or the 



pulsar wind shock or the shock in the jet or at the jet base. 
In the case of Cyg X-3 binary system, the jet scenario is 
more likely due to clear evidences of the existence of jets in 
this source (Marti et al. 2001, Miller- Jones et al. 2004). 

Let us consider the conditions in the jet which allow ac- 
celeration of electrons to a few TeV. The maximum energies 
of particles are determined by the balance between the en- 
ergy gain from the acceleration mechanism and energy losses 
on the radiation processes. In the case of the shock accel- 
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eration scenario, the energy gain rate can be parametrized 
by, 

Ace = £cE/R L « 10 13 £B eV s" 1 , (1) 

where E is the energy of electron, £ is the acceleration effi- 
ciency, -Rl = E/eB is its Larmor radius, B is the magnetic 
field inside the acceleration region (in Gauss), e is the elec- 
tron charge, and c is the velocity of light. The synchrotron 
energy loss rate of these electrons during acceleration pro- 
cess is given by, 

P syn ( 7 ) = 4a T c(7 B 7 2 /3 » 7.4 x 10 _4 B 2 7 2 eV s" 1 , (2) 

where (7b = B 2 /8iv is the energy density of the magnetic 
field, 7 is the Lorentz factor of electrons, and ctt is the 
Thomson cross section. By comparing Eq. \T\ and Eq. [21 we 
obtain the limit on the maximum energies of electrons ac- 
celerated at the shock, 

ESL « 60(£/B) 1/2 TeV. (3) 

Electrons can reach TeV energies (e.g. > 3 TeV), provided 
that the magnetic field in the acceleration site fulfills the 
condition (obtained by reversing Eq. [3} , 

B < 400£ G. (4) 

On the other hand, acceleration process of electrons can be 
also saturated by their IC energy losses in the strong ra- 
diation from the stellar companion (i.e. the hot WR type 
star). The IC energy losses in the Thomson (T) regime in- 
creases fast with energy. However, in the Klein Nishina (KN) 
regime, they increase slowly, depending only logarithmically 
with energy. The critical electron Lorentz factor, at the bor- 
der these two regimes is, 

1t/kn = mc 2 /3k B T ps2x 10 4 /T 5 , (5) 

where m is the electron rest mass, &b is the Boltzmann con- 
stant, T = 10 5 T5 K is the surface temperature of the com- 
panion star. Unless the acceleration process is saturated in 
the T regime, it is also not likely saturated in the KN regime 
since the energy losses in the KN regime depend weakly 
(only logarithmically) with electron energy. We evaluate the 
conditions under which acceleration of electrons can be sat- 
urated by IC energy losses in the T regime. The energy loss 
rate on this process is given by 

Ac(7) = 46T T cC/ rad 7 2 /3 « 1.3 x 10 4 T 5 4 7 2 /r 2 eV s"\ (6) 

where the energy density of radiation is I7 r ad = 4.7 x 
10 17 T 4 /r 2 eV cm -3 , and r is the distance from the com- 
panion star in units of stellar radius. The IC energy losses 
at the border between T and KN regimes can be approxi- 
mately estimated by putting Jt/kn into the above formula, 

P?c(7t/kn) ~5 x 10 12 T 5 2 /r 2 eV s" 1 . (7) 

Electrons can reach energies above mc 2 7T/M, when their 
energy gain rate is clearly above the IC energy loss rate at 
the border between the T and KN regimes. By comparing 
Eq. [T] and Eq. [7] we obtain the lower limit on the magnetic 
field strength at the acceleration region for which electrons 
can reach TeV energies, 

B > r 5 a /(2fr 2 ). (8) 



This limit is ~ 1/(10£) G for the case of the Cyg X-3 bi- 
nary system in which the temperature of the companion 
star is estimated on ~ 10 5 K and the distance of the com- 
pact object from the star is r — 2.25. The above limits on 
the magnetic field strength in the acceleration region which 
allow acceleration of electrons to 3 TeV, can be converted 
to the lower limit on the acceleration efficiency, £ > 0.016 
(from comparison of Eq. [4] and Eq.[8]). We conclude that in 
the case of the shock acceleration scenario (operating within 
the jet), electrons can be in principle accelerated to a few 
TeV, provided that the acceleration mechanism is efficient 
enough. Note, that £ in the Crab Nebula is expected to be 
not far from unity and the theoretical estimates give the val- 
ues of £ close to ~ 0.1 in the case of relativistic shocks (e.g. 
Malkov & Drury 2001). Applying the limiting case £ ~ 1, 
we estimate the magnetic field in the acceleration region of 
the jet on Bj ~ 0.1 — 400 G (able to accelerate electrons to 
a few TeV). Such magnetic fields in the jet of microquasar 
seem to be characteristic at distances from its base compa- 
rable to the dimensions of the companion star, i.e. acceler- 
ation process still occurs inside the binary system (see e.g. 
Bednarek 2006a, Bosch-Ramon, Romero & Paredes 2006). 
The high energy electrons and 7 -rays can be also injected 
within the binary system in other acceleration mechanism, 
e.g. as a result of electron acceleration in the reconnection 
regions in which case synchrotron losses may become negli- 
gible (e.g. Schopper, Lesch & Birk 1998), or as a secondary 
products of hadronic interactions occurring between rela- 
tivistic hadrons and the background matter inside the bi- 
nary system or at the vicinity of the compact object (e.g. 
Romero et al. 2003). 

TeV electrons produce first generation of 7-rays in the 
IC process by scattering anisotropic radiation from a com- 
panion star. These 7-rays are absorbed in the interaction 
with this same radiation field. As a result, 7-ray spectrum 
is formed in the IC e pair cascade process occurring in 
the whole volume of the binary system. The efficiency of 
the cascade process (the final 7-ray spectra escaping to- 
wards the observer) strongly depends on the location of 
the observer in respect to the source of primary electrons 
and the companion star. A basic prediction of such model is 
the anticorrelation of the GeV and TeV 7-ray emission with 
the orbital period of the binary system (see Bednarek 2000, 
Bednarek 2006b) . Such features have been recently observed 
in the GeV-TeV 7-rays from LS 5039 and LSI +61 303 (Abdo 
et al. 2009b,c). 

We apply this successful general cascade scenario to the 
binary system Cyg X-3, in order to investigate its likely 
behaviour at sub- TeV 7-ray energies. Such predictions can 
be very useful for planning future observations of Cyg X- 
3 with the Cherenkov telescopes. In order to control the 
number of free parameters of the model, we rather prefer 
to discuss general picture (without specifying any concrete 
model for particle acceleration). Therefore, we assume that 
somewhere inside the binary system a point like isotropic 
source of relativistic electrons appears. The location of the 
source is described by its distance from the companion star, 
JLy (measured in the plane of the binary system), and the 
distance, D (above the plane of the binary system). The 
source is located at the phase, <p (phase zero corresponds 
to the source behind the companion star in respect to the 
observer), and the observer is located at the inclination an- 
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gle i. The source accelerates electrons with the power law 
spectrum extending up to TeV energies. In our calculations, 
we apply the single power law differential spectrum for elec- 
trons, dN/dE cc E~ ZA , with the cut-off at 3 TeV. The spec- 
tral index has been chosen in order to be consistent with 
the differential spectral index of the GeV 7-ray emission 
observed by the Fermi satellite (Abdo et al. 2009a). With 
these assumptions on the source of relativistic electrons, we 
analyse the "pure" IC e ± pair cascade in the radiation of a 
massive star characterised by the parameters of the WR star 
observed in Cyg X-3. We perform cascade calculations un- 
der specific assumptions. At first, the local isotropization of 
secondary cascade pairs has been assumed. Such assump- 
tion is valid provided that the Larmor radius of electrons is 
much smaller than the mean free path, Aic, for their IC 
energy losses (Bednarek 1997), i.e. 

2ttR l < Aic. (9) 

The above condition is the most restrictive for electrons with 
energies at the border between the T and KN regimes (in 
the KN regime Rl/\™ « const). We estimate A^/ KN = 
cE/Piq^t/kn) ~ 6 x 10 7 r 2 /if cm, for electrons with en- 
ergy -Et/kn — m 7T/KN- Then, the condition given by Eq. 9 
can be transferred to the limit on the magnetic field in the 
region of cascade development, 

Bcas > 3.6TI /r 2 G. (10) 

For other energies of electrons than E T / KN , the limit is less 
restrictive. 

It is assumed that primary electrons and secondary lep- 
tons cool completely in the IC process up to 500 MeV. 
Note that more realistic scenario (which include a num- 
ber of possible additional processes) can be very compli- 
cated. It will need many free parameters for correct descrip- 
tion. Unfortunately, such parameters are at present not well 
known preventing detailed analysis of the complite cascade 
picture at this stage. For example, we neglected possible 
effects due to the strong magnetic field in the region of 
the cascade development such as synchrotron energy losses 
and partial re-distribution of directions of leptons by the 
ordered magnetic field (particles following local magnetic 
field). The importance of these effects were already consid- 
ered in (Bednarek 1997, see Eq. 3 and Fig. 4). It is possible 
to estimate the magnetic field strength in the cascade re- 
gion (i.e. the region around the companion star) for which 
the synchrotron energy losses of secondary cascade e ± pairs 
can be neglected in respect to IC energy losses. In the T 
regime, it is enough to compare the energy density of the 
magnetic field with the energy density of stellar radiation 
(Us = Kad). Then, the upper limit on the magnetic field 
strength in the cascade region is, 

B c T as < 4.3 x 10 3 T|/r G. (11) 

However, the limit on the magnetic field is much more re- 
strictive for electrons scattering radiation in the KN regime. 
By analysing dependence of the energy losses of electrons on 
these two radiation processes, we conclude, that with good 
approximation the limit on B in the KN regime should scale 
as B™ ~ BJssIEt /kn / E) ■ Therefore, for e pairs with en- 
ergies of 3 TeV, the limit on the magnetic field strength in 
the cascade region is of the order of ~ 10 G. Although the 
magnetic field around the companion star drops fast from 



its surface (close to the surface in the dipole region as oc r~ 3 
and in the region dominated by the wind as oc r -2 ), it may 
happen that synchrotron process may extract some addi- 
tional energy from the cascade e ± pairs. However, this will 
not influence the final conclusions of our paper on the sub- 
TeV 7-ray emission escaping from the binary system. The 
synchrotron losses will result in even lower level of the 7-ray 
fluxes than predicted by our "pure" IC e ± pair cascade sce- 
nario. Note, that synchrotron energy losses have been taken 
into account in some versions of the IC e pair cascade sce- 
nario (e.g. Sierpowska & Bednarek 2005, Bednarek & Gio- 
vannelli 2007 and also Bosch Ramon et al. 2008). 

We also neglected the adiabatic energy losses of leptons 
and non-local production of secondary 7-rays by secondary 
leptons (i.e. efficient escape of leptons from their place of 
origin). These two processes are difficult to treat reason- 
ably with a limited number of free parameters. We note that 
all the effects mentioned above can be considered in detail 
when precise multiwavelength results are available allowing 
realistic constraints of the free parameters describing these 
processes in the case of specific binary system. 

In the following subsections we calculate the 7-ray light 
curves at energies above 1 GeV, 100 GeV, and 330 GeV. 
We show the 7-ray spectra at specific phases of the injec- 
tion source. The results are discussed for the circular and 
eccentric orbits of the source of relativistic electrons. 

3.1 Gamma-ray light curves 

We calculate the cascade 7-ray light curves expected for two 
inclination angles (i = 30° and 60°) of the binary system 
Cyg X-3 assuming the point-like isotropic source of rela- 
tivistic electrons. Electrons have the power law spectrum as 
described above. The fluxes of 7-rays are obtained above 
three energies (1 GeV, 100 GeV, and 330 GeV), character- 
istic for the satellite observations and Cherenkov telescopes. 
In the first set of calculations (shown in Fig. 3), the injec- 
tion source is at a circular orbit but at different distances 
from the companion star. As expected, the GeV 7-ray light 
curves show clear modulation with the orbital period of the 
binary system with the maximum at phases when the injec- 
tion source is behind the companion star. For small inclina- 
tion angles (i = 30°), the 7-ray light curves are very similar. 
Note that the angle between the injection source and the ob- 
server does not change with the phase of the binary system 
for different radii of the orbit. This almost independence of 
GeV 7-ray light curve on the distance from the companion 
star is a consequence of a type of the IC e ± pair cascade 
considered in this paper (see Sect, above). Note that elec- 
trons interact mainly with soft photons coming from the 
direction of the companion star. When electrons are located 
farther from the star then the radiation field becomes more 
anisotropic and the first generation of 7-rays is more focused 
towards the companion star. These cascade 7-rays interact 
with the soft photons on average at similar distances from 
the companion star. As a result, most of the GeV cascade 
emission comes from the vicinity of the companion star. This 
effect is also responsible for the so called focussing effect of 
7-rays by a strong radiation of the companion star noted 
in Bednarek (2000). 7-ray modulation pattern changes with 
distance from the companion star for large inclination an- 
gles (i.e. 60°, compare the two upper figures in Fig. 3). Dif- 
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ferences are clearly seen for the phases when the injection 
source is in front of the companion star (phases close to 0.5). 
This is a consequence of the inefficiency of the cascade in the 
outward direction from the companion star. Note that the 
injection source is hidden behind the companion star in the 
case of an orbit with the radius of 1.5-R* and the phases 
close to zero (see dashed curves in the upper Fig. 3). That's 
why, 7-ray flux drops drastically. But, still some 7-ray pho- 
tons emerge from the opposite site of the star due to the 
focussing effect mentioned above. 

The sub-TeV 7-ray light curves show clear anticorrela- 
tion with the GeV 7-ray light curves. They peak at phases 
when the injection source is in front of the companion star. 
Sub-TeV 7-ray light curves show strong dependence on the 
radius of the orbit since on average this emission originates 
in the cascade process occurring at larger distances from the 
companion star. Moreover, the range of phases, at which 
sub-TeV emission is significant, decrease with the radius of 
the orbit of the electron injection source. For example, at 
energies above 100 GeV, the emission is concentrated in 
the range of phases between ~ 0.3 — 0.7 for the distance 
Rj = 20.R* but only in the range 0.4 — 0.6, for the distance 
R-t = 2.25r*. We conclude that interesting fluxes of sub-TeV 
7-rays can only have a chance to emerge from the binary sys- 
tem for the phases when the optical depths are below unity 
(compare Fig. 3 with Fig. 1). Stronger sub-TeV emission is 
expected for larger inclination angles of the binary system. 

The dependence of the cascade 7-ray light curves on 
the distance from the plane of the binary system, for the 
radius of the orbit of the injection source of electrons equal 
to R 1 — 2.25-R* (also the distance of the compact object) is 
shown in Fig. 4. Note that the level of modulation of the GeV 
7-ray flux with the period of the binary system decreases 
with the distance from the orbital plane being negligible for 
distances > 10i?*. For the inclination angle of the binary 
system equal to i = 30°, also the sub-TeV 7-ray flux shows 
weak modulation at large distances above the plane of the 
binary system. The level of modulation of the GeV-TeV 7- 
ray signal increases significantly for larger inclination angles 
(see the case for i = 60°). As in the previous cases, the 
GeV and sub-TeV emission is clearly anticorrelated. Again, 
for compact orbits of the injection source, the sub-TeV 7- 
ray flux is limitted to the range of phases when the injection 
source is in front of the companion star. This range of phases 
increases with the distance above the plane of the binary 
system. 

In summary, our cascade calculations show that the 
GeV and sub-TeV 7-ray emission is expected to be strongly 
modulated with the period of the binary system for the in- 
jection source of electrons within a few stellar radii from the 
companion star. The GeV and sub-TeV 7-ray light curves 
show clear anticorrelation, with the maximum of GeV emis- 
sion when the injection source is behind the companion star. 
The level of modulation of the GeV and sub-TeV 7-ray fluxes 
significantly drops with the distance from the companion 
star. We conclude that detection of the strong modulation 
of the GeV 7-ray flux indicates that 7-ray emission comes 
from within the binary system. In such a case, the sub-TeV 
7-ray emission is expected on a rather low level, strongly 
modulated with the period of the binary system, but with 
the maximum flux anticorrelated with the GeV 7-ray flux. 



3.2 Gamma-ray spectra 

We also show detailed 7-ray spectra at different phases of 
the binary system for the geometrical cases discussed above. 
The spectra are normalized to the peak emission detected 
by the Fermi satellite equal to 2 x 10~ 6 ph cm~ 2 s _1 (Abdo 
et al. 2009a). For comparison, we show the level of 1% of 
the Crab Nebula emission (marked by 1% C.U.) at energies 
above 100 GeV. It is clear that in all of the cases shown 
in Fig. 5 (dependence on the distance from the compan- 
ion star) and Fig. 6 (dependence on the distance from the 
plane of the binary system), the sub-TeV 7-ray emission is 
clearly below the sensitivity level of the present Cherenkov 
telescopes (e.g. the MAGIC or the VERITAS). Only the 
telescope array with the sensitivity on the level of ~ 0.1% 
C.U. (above 100 GeV) will have a chance to detect the 7-ray 
signal from Cyg X-3. On the other hand, Cyg X-3 might be 
detected by the on-ground telescope with sensitivity of the 
order of ~ 1% C.U. provided that its energy threshold is as 
low as ~20-30 GeV. It is not clear at present whether the 
next stage of the MAGIC and HESS arrays will obtain such 
sensitivity limit. We conclude that detection of the sub-TeV 
7-ray signal from the Cyg X-3 binary (WR type compan- 
ion star within compact binary system) by the Cherenkov 
telescopes has to probably wait for the construction of the 
next generation Cherenkov telescopes such as the Cherenkov 
Telescope Array (CTA). 

3.3 Circular versus eccentric orbit in Cyg X-3 

The shape of the GeV 7-ray light curve measured by the 
Fermi satellite shows a modulation with the period of the 
binary system but also clear asymmetry. The increasing part 
of the 7-ray light curve takes clearly longer than the decreas- 
ing part. Also the asymmetry in respect to the phase 0.5 is 
visible (the phase zero is counted in the Fermi plot from the 
location of the compact object behind the companion star, 
Abdo et al. 2009a). Such asymmetry in the GeV 7-ray light 
curve may be expected in the case of an eccentric orbit of 
the compact object (injection source of relativistic particles) 
along the companion star. 

As an example, we perform calculations of the 7-ray 
light curves for different inclination angles of the binary sys- 
tem and specific parameters of the eccentric orbit of the 
electron injection source located at or close to the compact 
object: the angle of the periastron is equal to uj = 60° (mea- 
sured in respect to phase zero), the eccentricity e = 0.3 and 
the semimajor axis a = 2.257?*. In our modelling, we as- 
sumed that: (1) electrons are injected inside a point like, 
isotropic source and interact only with the anisotropic radi- 
ation of the companion star; (2) the injection rate of elec- 
trons does not depend on the distance between the injection 
source and the companion star. The expected 7-ray light 
curves above three different energies are shown in Fig. 7. 
Interestingly, for the above parameters, we obtain the GeV 
7-ray light curves which general shapes strongly resamble 
the 7-ray light curve measured by the Fermi. Note, that the 
level of the modulation of the 7-ray flux with the orbital 
period increases with the inclination angle, changing by a 
factor of ~5 for the inclination i = 30° and up to a factor 
of ~100 for the inclination angle i = 60°. At present, it is 
difficult to decide which inclination angle describes the ob- 
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Figure 3. The cascade 7-ray light curves expected for a point like isotropic source of relativistic electrons located on a circular orbit 
around the companion star in the plane of Cyg X-3 binary system for energies: E 1 > 1 GeV (upper panel), > 100 GeV (middle), and 
> 300 GeV (bottom). The phase zero corresponds to the location of the injection source of electrons behind the massive star. Electrons 
are injected with the differential power law spectrum, oc E~' iA , up to 3 TeV. The injection rate of electrons is independent on the phase 
of the binary system. Figures on the left are obtained for the inclination of the binary system equal to i = 30° and on the right for 
i = 60°. The radius of the orbit around the companion star is equal to: Ry = 1.5R+ (dashed curves), 2.25-R* (solid), 5-R* (dot-dashed), 
10-R* (dotted), 20-R* (triplc-dot-dashed). The companion star has the radius R+ = 2 X 10 11 cm and the surface temperature T* = 9 X 10 4 
K as expected for the case of Cyg X-3 binary system. 



servations better due to the lack of knowledge on the precise 
location of the baseline emission in the Fermi light curve. 

We also show the expected 7-ray light curves at sub- 
TeV energies for three inclination angles of the binary sys- 
tem. A clear unticorrelation of the GeV and sub- TeV emis- 
sion is also expected in the case of Cyg X-3 with the max- 
imum emission at phases ~0.3-0.4. The largest sub- TeV 7- 
ray fluxes are expected for large inclination angles. We also 
show the 7-ray spectra for specific phases of the binary sys- 
tem and compare them with the sensitivities of the present 
Chercnkov telescopes (around ~ 1% Crab Unit). The spec- 
tra have been normalized to the GeV peak emission reported 
by the Fermi. A clear break in the 7-ray spectra is observed 
between 10-100 GeV in most of these spectra. Its precise lo- 
cation depends on the phase of the binary system. However, 
these spectra are not able to reach the level of the 1% C.U. 
So then, they are undetectable by the present Cherenkov 
telescopes. 

Note that these calculations have been performed under 
the assumption that the injection source is in the plane of the 
binary system (identified with the compact object inside the 
binary). However, as we have shown above (for the case of a 



circular orbit), the source of relativistic electrons can not be 
located far away from the plane of the binary system since 
in such a case the modulation of the GeV 7-ray signal drops 
significantly being in contradiction with the measurements 
by the Ferrni-LAT telescope. 



4 DISCUSSION AND CONCLUSION 

We have applied the anisotropic IC e pair cascade model, 
developed for the 7-ray production in the TeV 7-ray bina- 
ries, to the X-ray binary system Cyg X-3, recently discovered 
as a source of variable and modulated GeV 7-ray emission. 
The 7-ray light curves and spectra are calculated for dif- 
ferent locations of the injection source of relativistic elec- 
trons inside the binary system following two basic models: 
(1) injection of 7-rays mainly within (or close to) the plane 
of the binary system as expected in the accreting neutron 
star model (Bednarek 2009); (2) injection of relativistic elec- 
trons from the place above the plane of the binary system 
as expected e.g., in the microquasar model in which par- 
ticles are accelerated in the jet launched from a compact 
object (e.g. Levinson & Blandford 1996, Georganopoulos et 



On the possibility of sub-TeV Gamma-ray emission from Cyg X-3 9 



> 1 GeV 




dj O 



0.2 0.4 




>100 GeV 



>5l 





> 1 OO GeV 














y 

~ ✓* 









0.4- 0.6 O. 

>330 GeV 



0.4 0.6 

>330 GeV 



0.4- 0.6 
phase 



0.4- 0.6 
phase 



Figure 4. As in Fig. 3 but for the injection source of electrons located at the distance of the compact object (a = 2.25-R*) but above 
the plane of the binary system: D = 1_R* (solid), 3-R* (dot-dashed), 5R* (dotted), 10-R* (triplc-dot-dashcd), and 20-R* (dashed). 



al. 2002, Romero et al. 2002). At first, calculations are per- 
formed for circular orbits of the injection source in order 
to investigate their basic features. For reasonable parame- 
ters of the Cyg X-3 binary system, we obtained the GeV 7- 
ray light curves which have general features consistent with 
the observations of this source by the Fermi-LAT telescope. 
Moreover, it is found that the sub-TeV 7-ray emission is 
not expected to be produced on the level allowing its de- 
tectability by the present Cherenkov telescopes. This con- 
clusion has been confirmed by some example calculations 
for a more complite scenario in which the injection source 
of electrons is located on an eccentric orbit around the com- 
panion star. The lack of detectable sub-TeV 7-ray emission 
from Cyg X-3 is a consequence of huge optical depths for 
7-ray photons propagating within the binary system. Such 
large optical depths are the consequence of a small compact- 
ness of the binary system and high surface temperature of 
the companion star (~ 10 5 K). Detection of the 7-ray signal 
from sources similar to Cyg X-3 with on-ground telescopes 
will be only possible with the next generation of Cherenkov 
telescopes which are expected to have sensitivity on the level 
of ~ 0.1% of Crab Units or by the telescopes which will be 
able to perform observations with the sensitivity of ~ 1% 
C.U. in the 20-30 GeV energy range. Such sensitivities are 
planned for the future Cherenkov Telescope Array (CTA). 
Note however, that in the case of Cyg X-3 another complica- 
tion is introduced by the transient nature of the GeV 7-ray 
emission (a few day outbursts, Tavani 2009). This will ad- 



ditionally lower the chances for detection of Cyg X-3 in the 
TeV energies even with the CTA telescopes. The chances for 
detection can significantly rise if the 7-ray spectrum in the 
GeV energies is natter during some outbursts than reported 
recently (differential spectral index -2.7). Therefore, investi- 
gation of a larger number of GeV 7-ray outbursts from Cyg 
X-3 by the AGILE and the Fermi telescopes will be of great 
importance for planning future sub- TeV 7-ray observations 
with the Cherenkov telescopes. 

In the case of the microquasar model for the 7-ray pro- 
duction in binary systems two jets, propagating above and 
below the plane of the binary system, are expected. The 
calculations shown here concern only the 7-ray production 
from one of the jets (i.e. the jet propagating into the hemi- 
sphere containing the observer). For large inclination angles 
of the binary system, the contribution of the jet and counter- 
jet becomes comparable but for small inclination angles the 
counter-jet contribute mainly to the GeV energy range due 
to large angle between the injection source of electrons and 
the observer (measured from the companion star) . As a re- 
sult, normalization to the Fermi spectrum at the GeV en- 
ergies to the sum of the 7-ray emission from the jet and the 
counter-jet in the case of small inclination angles of the Cyg 
X-3 binary system additionally decreases estimated in this 
paper 7-ray emission at sub-TeV energies. 

Finally, we would like to comment of the possible in- 
fluence of the assumptions in considered model on the pre- 
dicted 7-ray light curves and spectra from Cyg X-3. 
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Figure 5. The 7-ray spectra produced by electrons injected from a point-like isotropic source in the orbital plane of the binary system 
but at different distances from the companion star: i? 7 = 1.5.R* (a and f), 2.25-R* (b,g), 5i?* (c,h), 10i?* (d,i), and 20R* (e,j). Primary 
electrons have a power law spectrum with differential spectral index —3.4 extending up to 3 TeV. Specific curves show the spectra for 
different phases of the injection source: ip = 0. (dashed curves), O.f (dot-dashed), 0.3 (dotted), 0.4 (triple-dot-dashcd), and 0.5 (solid). 
The phase zero corresponds to the location of the injection place behind the massive star. Figures on the left arc obtained for the 
inclination of the binary system equal to i = 30° and on the right for i = 60° . The thin solid line shows the GeV peak emission observed 
by the Fermi and the thin dotted line marks the sensitivity of the modern Cherenkov telescopes which is on the level of 1% Crab Unit. 
The companion star has the radius = 2 X 10 11 cm and the surface temperature T* = 9 X 10 4 K, as expected for the case of Cyg X-3 
binary system. 
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Figure 6. As in Fig. 5 but for the parameters as in Fig. 4. 



(i) Synchrotron energy losses of primary and secondary 
leptons: This process will additionally extract energy from 
leptons in respect to considered pure IC cascade process (e.g. 
Bosch Ramon et al. 2008). It is expected that with impor- 
tant synchrotron energy losses the calculated 7-ray fluxes 
specially at larger energies > 100 GeV becomes even lower 
due to the dominance of the synchrotron energy losses of 
electrons over their IC energy losses in the Klein-Nishina 
regime. 

(ii) Incomplete cooling of leptons and their escape from 
the production place: First process will result in deficiency 



of 7-ray flux at lower energies (GeV range). In order to be 
consistent with the spectral index of the 7-ray spectrum 
even a steeper spectrum of primary electrons should be ap- 
plied. However, this will result in lowering also the 7-ray 
flux above ~ 100 GeV. The influence of the second process 
is difficult to estimate since it is even not clear in which di- 
rection secondary leptons can be driven from the their site 
of origin. The effects of advection in the fast stellar wind 
seems to be negligible due to relatively large advection time, 
Tadv ~ Ry/vw (where t) w is the stellar wind velocity which 
in the case of WR stars is of the order of ~ (2 — 3) x 10 km 
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Figure 7. The example 7-ray light curves for the eccentric orbit 
of the injection source of electrons with the angle of the periastron 
passage ui = 60°, the eccentricity e = 0.3, and the scmimajor 
axis a = 2.25-R*. At phase zero, the injection source (a compact 
object) is behind the companion star. 7-rays are produced in the 
IC pair cascade process initiated by primary electrons (with 
the differential power law spectrum <x E~ 3,4 extending up to 
3 TeV) injected from a point like source. The light curves are 
shown for 7-rays with energies above 1 GeV (top panel), 100 
GeV (middle), and 330 GeV (bottom). The inclination angle of 
the binary system is equal to: i = 30° (dashed curve), 45° (solid), 
60° (dotted). 



s _1 ), in respect to the cooling time of electrons on the IC 
process. The IC cooling time is t^ c Kt 40r 2 /(T57) s, where 
the dilution factor of radiation from the stellar surface can 
be approximated by r = Ry/R*. Even for electrons with en- 
ergies 500 MeV, injected at the distance, R-y = 10.R*, from 
the star with surface temperature T = 10 5 K, r a( jv ~ 10 4 
s and Tjc ~ 4 s. So then, tic << i~adv even for electrons 
scattering soft radiation in the Klein-Nishina regime with 
the largest considered energies equal to 3 TeV. 

(iii) Izotropization of cascade leptons versus their propa- 
gation along the local magnetic field lines: The linear cascade 
through the radiation field of the companion star has been 
also considered as the limitting case when considering the 



importance of the IC e pair cascade contribution to the 7- 
ray spectrum escaping from the binary system (e.g. Aharo- 
nian et al. 2006, Cerutti, Dubus & Henri 2009). We discussed 
here scenario in which electrons are locally isotropized. The 
general conditions for local izotropization of cascade leptons 
has been discussed in Sect. 3. The isotropization process 
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Figure 8. The 7-ray spectra produced at different phases of the 
injection source of rclativistic electrons in the binary system Cyg 
X-3. The parameters of the orbit and the spectrum of primary 
electrons arc as in Fig. 7. The inclination angle of the binary 
system is fixed on: i = 30° (a), 45° (b), 60° (c). The spectra for 
phases ip = 0. (triple-dot-dashed curve), 0.2 (dashed), 0.4 (solid), 
0.6 (dot-dashed), 0.8 (dotted). 



occurs efficiently (based on the comparison of the Larmor 
radius with the mean free path for IC process) even in a 
relatively low magnetic field. However, exact propagation of 
leptons in the ordered component of the magnetic field will 
depend on their injection angles in respect to the magnetic 
field direction. Magnetic focussing may result in production 
of strong enhanced 7-ray emission at specific directions on 
the sky directly unrelated to directions of the first genera- 
tion of cascade 7-rays. These effects have been studied in 
detail by Sierpowska & Bednarek (2005). 

(iv) Phase dependent injection rate of electrons: The in- 
fluence of this effect is also very difficult to consider without 
detailed knowledge on the mechanism of acceleration process 
and generation of energy. It might be expected that both 
these processes could depend on the distance between the 
compact object and the companion star. However, whether 
the higher accretion rate (or stronger stellar wind) will in- 
crease or decrease power converted to relativistic electrons or 
how this will influence their maximum energies is at present 
not clear. Therefore, the only reasonable approach seems to 
keep these basic processes as independent on the distance 
between the stars (i.e. independent on the phase of the bi- 
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nary system). These effects can be studied when detailed 
multi-wavelength observations of 7-ray binary systems be- 
come available. 

(v) Possible absorption of 'y-rays in the X-rays: We pos- 
tulate that 7-rays, produced in terms of the jet model, orig- 
inate at distances from its base which are comparable to the 
dimensions of the binary system. The radiation from the ac- 
cretion disk and the inner part of the jet can produce strong 
target even for hadrons (see e.g. models by Levinson & Wax- 
man 2001, Bednarek 2005, or more recently by Romero & 
Vila 2008). So then, also 7-rays produced in the inner region 
of the jet are strongly absorbed. However, 7-rays, produced 
at large distances from the base of the jet, suffer the radia- 
tion form the inner disk and/or jet significantly diluted. At 
distances of the order of thousands larger than the inner disk 
radius, this radiation field is diluted by a factor proportional 
to the square of the distance, i.e. by a factor of ~ 10 6 . We 
conclude that such strongly diluted radiation field from the 
inner disk and/or jet will not be able to efficiently absorb 
the high energy 7-rays produced farther away from the base 
of the jet. 
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